Aharonov-Casher oscillations of spin current through a multichannel mesoscopic ring 
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The Aharonov-Casher (AC) oscillations of spin current through a 2D ballistic ring in the pres- 
ence of Rashba spin-orbit interaction and external magnetic field has been calculated using the 
semiclassical path integral method. For classically chaotic trajectories the Fokker-Planck equation 
determining dynamics of the particle spin polarization has been derived. On the basis of this equa- 
tion an analytic expression for the spin conductance has been obtained taking into account a finite 
width of the ring arms carrying large number of conducting channels. It was shown that the finite 
width results in a broadening and damping of spin current AC oscillations. We found that an ex- 
ternal magnetic field leads to appearance of new nondiagonal components of the spin conductance, 
allowing thus by applying a rather weak magnetic field to change a direction of the transmitted spin 
current polarization. 

PACS numbers: 73.23.-b, 85.75.-d, ,75.76.-fj, 71.70.Ej 



I. INTRODUCTION 

The Aharonov-Casher effect is a spectacular demon- 
stration of the fundamental role of the spin-orbit inter- 
action (SOI) in electronic transport. This effect is a non- 
Abelian analog of the Aharonov-Bohm effect. An elec- 
tronic wave entering a two dimensional ring is splitted 
into two waves traveling trough the upper and the lower 
arms and interfering on the exit from the ring. Due to 
SOI the spinor components in each of the waves obtain 
the relative phase shift. This shift, in its turn, gives rise 
to a destructive or constructive interference pattern in 
the transmittance probability, that results in a number of 
oscillation effects on electron transport parameters. No- 
tably, that in semiconductor heterostructures SOI can be 
varied through the gate voltage manipulation, suggesting 
interesting opportunities for practical applications of this 
effect in spintronics. 

A simplest model to study the AC effect is a one- 
dimensional (ID) or one-channel ring with none or few 
scatterers of electrons. Aronov and Lyanda-Geller^ stud- 
ied oscillations of the magnetoconductance in a ID ring 
in the presence of Rashba^ SOI. Ya-Sha Yi et. al.^ con- 
sidered a joint effect of the Zeeman coupling, magnetic 
flux and SOI on the conductance of a ID ring beyond 
the adiabatic approximation employed in Ref. 1. Nitta 
et. al.'^ noted that SOI alone, without the magnetic field, 
can cause oscillations of the ring electric conductance. 
Meijer et al.^ pointed out that the Hamiltonian used ear- 
lier in Ref. 1 was not quite correct and, in particular, was 
not Hermitian. Using the correct Hamiltonian Frustaglia 
and Richter® revised the expression for the conductance 
found in Ref. 4. 

In disordered mesoscopic systems the AC effect is 
strongly modified. Multiple impurity scatterings lead to 
averaging out of some strong oscillations that were pre- 
sented in ideal ID rings. There is a similarity to AB 
effect where the fundamental h/e peak in the Fourier 
spectrum of magnetoconductance vanishes and is sub- 
stituted for h/2e Al'tshuler, Aronov, and Spivak^ oscil- 
lations. Mathur and Stone* have demonstrated that a 



similar weak localization effect takes place in the case of 
AC oscillations in a thin diffusive ring with Rashba SOI. 
They have shown that the average conductance period- 
ically varies as a function of the SOI coupling constant 
with a period that is 1/2 of the ideal ring conductance 
oscillations. This prediction have been confirmed exper- 
imentally by Koga et al.^ who observed that the ampli- 
tude oi h/2e magnetoconductance oscillations (AAS type 
oscillations) oscillates itself with varying the applied gate 
voltage and, correspondingly, the SOI strength. On the 
other hand, the fundamental AB, as well as AC oscil- 
lations show up in mesoscopic conductance fluctuations. 
The conductance correlator has been calculated in Ref. 10 
for a diffusive ring in the presence of Zeeman coupling 
and Rashba SOI. It was found that the amplitude of h/e 
AB peak shows oscillations with varying SOI, similar to 
those that have been observed for the h/2e peak in Ref. 9. 

Although ID, as well as diffusive models are useful 
to elucidate fundamental physical effects associated with 
SOI, they can not be fully applied to realistic semicon- 
ductor systems used in experiments. For example, an 
attempt to interpret the observed^^ splitting of the AB 
power spectrum within the diffusion model^"'^^ was not 
successful. In Ref. 11, as well as in other experimental 
works^'^^, the mesoscopic loops carry many channels. At 
the same time their sizes are comparable to the electron 
mean free path. Therefore, neither of the above theo- 
retical models can be applied. In this situation an ap- 
proach based on path integrals along classic trajectories 
can be fruitful. Such a method has been applied to trans- 
port in mesoscopic systems in a number of works^^^^^. It 
was also employed for calculation of the spin conduc- 
tance through a classically chaotic and regular cavities 
and rings with Rashba SOI^*'^^. In Ref. 18 it was done 
analytically by applying the method of trajectory aver- 
aging, while numerical simulations of path integrals were 
performed in Ref.l9. In both cases pronounced AC os- 
cillations of the spin conductance with varying SOI con- 
stant had been found. Their important distinction from 
oscillations of the electric conductance is that they ap- 
pear in the main semiclassical approximation, not involv- 
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ing weak localization or other quantum corrections. 



II. THE MODEL 



While in a multichannel loop studied in Ref. 18 the 

electron motion was two-dimensional, the AC phase ac- 
cumulation had an effectively one-dimensional character. 
The reason is that in a thin enough loop SOI causes only 
a small variation of spinor amplitudes during particle mo- 
tion along any straight segment of a trajectory. In the 
leading approximation ignoring AC phase fluctuations as- 
sociated with a finiteness of a loop width, a phase evo- 
lution depends only on a coordinate along the loop and 
finite width effects vanish. On the other hand, as follows 
from the Monte Carlo analysis^^, when a particle lifetime 
within the ring is long enough, the finite width effect on 
the amplitude and shape of AC oscillations is strong. 

In order to elucidate this problem we will calculate the 
spin conductance taking into account the finite width ef- 
fect in a classically chaotic multi-channel ring with the 
Rashba spin-orbit interaction and a uniform magnetic 
field applied perpendicular to the ring. A chaotic mo- 
tion of particles can be provided by non ideal boundaries 
of the ring, as well as by random potential variations 
inside it. Starting from an analog of the Landauer for- 
mula derived for the spin conductance in Appendix A, 
we apply the path integral method to this conductance. 
Within this approximation spin dependent transmission 
amplitudes for each of the classical paths decompose into 
a product of a spin independent transmission amplitude 
and a matrix determining evolution of spinor components 
along this trajectory. The expression for the spin con- 
ductance that is quadratic in these amplitudes should 
be double summed over the trajectories afterwards. Ig- 
noring weak localization and other quantum corrections, 
only the terms diagonal with respect to the trajectories 
are then retained. As a result, the spin conductance 
takes a form of a 3 x 3 matrix averaged over the trajec- 
tories. Based on known statistical properties of chaotic 
trajectories we will derive Fokker-Planck equations de- 
scribing spin evolution of a particle moving through the 
ring and analytically calculate the average of spin con- 
ductance matrix components. The finite width effects 
will be analyzed that show up in an additional broad- 
ening and decreasing of AC oscillations. The magnetic 
field, in its turn, results in appearance of spin conduc- 
tance components that were equal to zero in the absence 
of the field, suggesting thus an opportunity to rotate the 
spin polarization on the exit from the ring. 

The outline of the paper is as follows. Section II con- 
tains a description of the model system we used in our 
theory. In section III an expression for the spin conduc- 
tance is obtained in the form of a sum over classical tra- 
jectories. In section IV the Fokker-Planck equation for 
the spin polarization distribution function is derived. On 
the basis of this equation the spin conductance averaged 
over chaotic trajectories is calculated. The discussion of 
results is presented in Section V. 



We consider spin transport through a 2D ring which 
is connected via two symmetrically placed leads to two 
reservoirs of electrons (see Fig. 1) and is subject to the 
magnetic field perpendicular to the ring plane. The latter 
gives rise to the Zeeman interaction 
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where e is the electron charge, H is the magnetic field, g 
is the g-factor, m* is the effective electron mass, and c is 
the light velocity. 

The Rashba spin-orbit interaction is assumed to take 
place only in the range of the ring. It has the form 



5Hi 



, (T • [p X Zj 



(2) 



Here ago is the SOI constant, the vector & consists of the 
Pauli matrices ax, a-y and iTz, p is the momentum oper- 
ator and z is a unit vector parallel to the Z-axis. We 
assume the hard wall reflection of electrons from the ring 
boundaries. Since the particle spin is conserved upon 
such a reflection, Eqs. (1) and (2) determine a spin dy- 
namics in the ring. 

The reservoirs are assumed to be in a local thermody- 
namic equilibrium with a given polarization (magnetiza- 
tion). For simplicity we assume that the left reservoir 
is polarized (along a unit vector u), while the right one 
is unpolarized. The polarization of the left reservoir is 
characterized by the chemical potentials Hl^ = I^l + ^ 
and = Ml — for spin- up and spin-down (relative 
to I/) electrons, respectively. For the right reservoir we 
assume, in its turn, that /z_r-|- = /x/jj, = Hl- This situation 
is of particular interest for our further analysis, because 
establishing of thermodynamic equilibrium between spin 
subsystems in the reservoirs will be accompanied by the 
spin current, while the electric current will be absent. In 
the linear response regime the spin current in the right 
lead can be expressed (see Appendix A) as 



(3) 
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FIG. 1: Geometry of the problem under consideration. The 
Z-Bxis points to the reader. W is the leads width, d is the 
ring arms width, and a is the radius of the ring. 
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where i and j take the values x,y,z, and uj is the j-th 
component of u. 

III. SEMICLASSICAL APPROXIMATION 

We assume that the leads connecting the ring with 
reservoirs are ideal conductors with a constant cross sec- 
tion. This siiggcsts the use of the Landaiicr approach 
for calculation of the spin current. The latter, however, 
has an important distinction from the electric current, 
because it does not conserve in the region with SO inter- 
action. Let us assume, for example, that such a region 
is ideally transparent. Even in this case the transmitted 
spin current will not be equal to the incident one. Alter- 
natively, we may consider a polarized reservoir connected 
via one lead to a region with SO interaction. In the sta- 
tionary case the electric current through this lead will be 
zero. At the same time , the spin current will be finite, 
because spin polarizations of incident and reflected elec- 
trons can be different. This sort of spin transport has 
been considered in Ref. 18 where a ring played the role 
of the region with SO interaction. 

So, we see that Landauer formula should be modified 
to take into account effects of nonconcerving spin. In 
general case it is done in Appendix A. For a particular 
set up considered in this work the expression for the spin 
conductance takes a simple form, somewhat similar to 
the Landauer formula 

" 4^ ^ {*P^p(/^fl) Vp(Mfl) } (4) 



pp' 



X 2 matrix composed of 
tt Ai f,,\ A1- 



Here tpip(iifi) is the 2 
the transmission amplitudes tplp{iJ,R),tpfp{fiR) , ^^/^(Mfi) 
tp^pil^n) from the channel p of the left lead to the chan- 
nel p' of the right lead. The summation is performed 
over open channels. Note, that this formula implies that 
for calculation of all spin conductance components it is 
not sufficient to calculate the spin-resolved transmission 
coefficients " = J2p'p Kp'p P (studied, for example, in 
Refs. 6,20). This is readily seen, for example, by inspect- 
ing the expression 
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pp 



which follows from (4). 

To find the transmission amplitudes in (4) we will fol- 
low Ref. 21 (see also Ref. 22) for the spinless transmis- 
sion amplitude at the Fermi energy Ep. In the presence 
of spin degrees of freedom it allows the evident general- 
ization: 

tp'p = -ihy/v^ j dydy'u*p,{y')up{y)G{y'a' ,ya\EF) ■ 

(5) 

Vpi (vp) is the longitudinal velocity in the p'-th (p-th) 
channel of the right (left) lead. G is the retarded 



Green's function. The integration is performed over 
cross-sections of the leads. We assume that the leads 
and ring arms are wide enough, so that there are many 
channels below the Fermi energy. This allows to apply 
the semiclassical approximation (see e.g. Ref. 14) to 
Eq. (5). Within this approximation the path integral 
expression for the Green's function is replaced by the 
sum of amplitudes corresponding to classical trajecto- 
ries traversing the ring (details of this procedure may be 
found in Ref. 17). Integration over y and y' is performed 
by the stationary phase method (large parameter in the 
exponent is the number of open channels in the leads). 
The result is 



(6) 



The label s enumerates the trajectories that enter the 
ring at the angle 6 = isin""'^ {pTT/kpW) relative to the 
a;— axis and exit it at 9' = ±sm~^ {p'Tr/krW), where 
W and kp are the leads width and Fermi wave number, 
respectively. to(s) is the spin independent transmission 
amplitude corresponding to the s— th classical trajectory. 
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sgn [Bs] sgn J A, 



X exp <^ ikpLs + ikp (sin^gi/s — sin^^y^) 



(7) 



Here Nm = int{kpW/Tr} is the number of open channels, 
Ls is the length of the s-th trajectory, yg {y'g) stands 
for the 2/-coordinate of the entrance (exit) point of s-th 
trajectory, A is the vector potential corresponding to the 
magnetic field applied to the ring. Other quantities in (7) 
are 
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dO{y,y') 



dy{e,e') 



dy 
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d9'{9,y') 



dy' 



where jis is the Maslov index^"''^'' and O is the Heavi- 
side step function. The matrix 5" " in (6) determines 
an evolution of the spin state along s-th trajectory. It 
should be noted that Eq. (6) has been derived assum- 
ing that classical trajectories do not depend on the spin 
dynamics. This allowed to write the terms entering into 
the sum in Eq. (6) in the form of a product of spin 
dependent and spin independent parts. In fact, this as- 
sumption means that in the leading semiclassical approx- 
imation we ignore a difference between Fermi velocities 
corresponding to spin-split subbands. It can be done, if 
during the time between two consecutive collisions with 
ring boundaries, a divergence of two wave packets belong- 
ing to these subbands will be much less than the electron 
wavelength. The corresponding condition can be written 
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as Lso = fi/oiso '77,* ^ d where Lso is the spin-orbit 
length that measures the SOI strength and d is the ring 
width. 

For a spin dependent Hamiltonian consisting of the 
two terms represented by Eqs. (2) and (1), the evolution 
operator can be expressed as 
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exp 
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dt ttsoPF CT ■ Tit 



hojH 



(8) 

Here T is the time ordering symbol, Pf is the Fermi mo- 
mentum, rtt is the unit vector parallel to p x z, and p 
(IpI ^ Pf) is the electron momentum. Note, that direc- 
tion Tit of the effective magnetic field generated by the 
SO interaction changes its sign when a particle reverses 
its motion direction. 

Substituting (6) into (4) we obtain 



9ji = ^X!X!*o(s)*o(u)Tr|ai5+aj5„| . 



(9) 



p\p s,u 



In this equation each semiclassical amplitude to{s) con- 
tains a phase factor exp {27rzXs/A}. Since the path 
lengths Lg and L„ of trajectories in Eq. (9) are much 
longer than the electron wavelength A, the terms with 
,s 7^ u oscillate rapidly even with a small variation of the 
particle energy, as well as slight change of the loop shape 
and/or impurity positions. In an experimental situation 
it can also be gate voltage variations and magnetic field 
switching^ ^. On the other hand, the diagonal terms with 
,s = u do not oscillate. If one is not interested in meso- 
scopic fluctuations of the spin conductance, or quantum 
corrections to it, only the terms with s = u have to be 
retained^^. On the basis of the crgodic hypothesis of 
Ref. 23 this procedure may also be treated as averaging 
over a random ensemble of the rings. 

Thus, from (9) the averaged spin conductance {gij) is 
obtained as 



(10) 



where 



K!j = ^Tv{aAajSl} . (11) 

The electrical conductance G after such averaging takes 
the form 



{G) = {e/nn)J2Ms) 



(12) 



where the factor 2 accounts for the spin degrees of free- 
dom. It should be noted that (G) given by (12) does not 
depend on the SO interaction strength. On the other 
hand, oscillations of G with the varying Rashba constant 
have been predicted in Ref. 1. This distinction can be 
explained by importance of quantum effects in an ideal 
ID ring considered in Ref. 1, while such effects are small 



in a disordered multichannel system that we study here. 
They can be taken into account as weak localization cor- 
rections to AC oscillations similar to those studied for 
diffusive rings. * 

Let us now consider the transparency 



p'p I 



for a spinless particle. Here tp>p is the transmission am- 
plitude from the channel p on the left to the channel p' 
on the right. So, the transparency N is normalized to the 
number of open channels N^. Applying configurational 
averaging to the above expression we obtain it as a sum 
over trajectories: 



(13) 



Hence, the ratio 



Pis) 



Nci 



can be identified with the probability that an electron 
chooses the s-th trajectory to pass the ring. Expressing 
|to(s)|^ from the last formula and substituting it into (10) 
we obtain 



(14) 



where {■)s denotes averaging over trajectories with the 
probability P{s) and 2e^5o = e^Nci/irh is a classical con- 
ductance of the ring^^ . 

We assume that the classical motion of an electron in- 
side the ring is chaotic. The chaotic dynamics is provided 
by small scale bumps and other irregularities on the ring 
boundaries, while macroscopically the ring preserves its 
regular shape. For chaotic trajectories after long enough 
time an electron "forgets" through which of the leads it 
entered the ring. Together with an assumption that the 
leads are symmetric this allows to conclude that probabil- 
ity for a particle to be reflected is equal to the probability 
to be transmitted. So, Nd = iVm/2 and go = Nm/Airh. 

Formula (14) will be used for calculation of the spin 
conductance (gij). However, another physically more 
transparent representation of (gij) can be suggested. As 
shown in Appendix B, K'^^is the i-th component of the 
electron polarization at the end of the s-th trajectory, 
provided that at its beginning the electron was polarized 
along the j-th axis. Taking this into account one can 
introduce the effective polarization vector P^f f{s\ej) = 
[gis)/go] {K^j,Kyj,Klj] of electrons at the end of the 

s-th trajectory. Here g{s) = \tQ{s)^ /2'Kh is the flux of 
spinless particles that pass the ring through the s-th tra- 
jectory. Eq. (10) then takes the form 



(15) 



5 



From comparison of (14) and (15) we see that {Kfj)s may 
be treated as i-th component of the effective polarization 
vector Peff{ej) on the exit from the ring, 



(16) 



Equation (15) means that the spin current on the exit 
from the ring is given by the sum of effective polariza- 
tions corresponding to each classical trajectory, times the 
flux of particles passing the ring. From this fact an im- 
mediate conclusion follows: if the resonance condition is 
satisfied, that is the effective polarizations corresponding 
to different trajectories on the exit from the ring are co- 
directional, then the spin current will have a maximum. 



IV. FOKKER-PLANCK EQUATION 

According to Eq. (14) calculation of the spin conduc- 
tance (gij) reduces to averaging of the polarization trans- 
formation matrix Kfj over trajectories. To perform this 
averaging the distribution function of K^^ is needed. One 
of the standard ways to find it is the following. We note 
that the quantities K^j in (11) correspond to the end of 
the s-th trajectory. It is evident however that one may 
extend definition (11) to any time instant t on the tra- 
jectory. Taking then a time derivative of (11) we arrive 
at the stochastic differential equation (of the Langevin 
type) for Kfj. This equation may be used to calculate 
drift and diffusion coefficients in the Fokker-Planck equa- 
tion. Solving then this Fokker-Planck equation one can 
find the desired distribution function. 



A. Dynamical equation for the spin S-matrix along 
a trajectory 

As a preparation step to our calculation we note that 
at H — eq. (8) for Sg can be written in the form of a 
contour integral over the s-th trajectory: 



5,(F = o) = r 



exp 



Lso 



dl ■ [z X &] 



(17) 



As any unitary operator acting in the spin space, Ss{H = 
0) can also be written in the form of the rotation operator 



^,(i/ = 0)=exp^z^| 



(18) 



where (Tjv is a projection of <t onto some unit vector 
N. An important consequence of (17) is that TV and 
^ depend only on the geometry of the s-th trajectory 
and don't depend on the dynamics of motion along the 
trajectory. This is the reason why 5' may be called a 
geometric phase^^. 

To extend the analysis of Ref. 18 and take into account 
finiteness of the ring width we divide the time interval 




FIG. 2: Replacement of the real trajectory with the fictitious 
one in the case of the ring of finite width. Both trajectories 
give rise to the same evolution operator Ss in the absence of 
the magnetic field. 



(0,i) into small subintervals Ati = ti — to = < 
ti < ■ ■ ■ < tn = t. After that S can be represented in the 
form (the index s is omitted when it doesn't lead to any 
confusion) 



S{t)^S{Atr,)S{Atn-i)...S{Ati). 



(19) 



The real trajectory is transformed next by adding to each 
i-th segment a path passed in direct and opposite direc- 
tions, as shown in Fig. 2. It is seen from (17), that paths 
passed twice in the opposite directions don't contribute 
to Ss{H = 0). Consequently, each term S{Ati) in (19) 
may be replaced without changing Ss with 



(20) 



Here b is the trajectory segment corresponding to the 
time interval Ati, o- and c connect the middle line of the 
ring with b, as explained in Fig. 2. Expanding exponents 
in (8) and (17) one obtains 



Sn 



1 



L 



rdl 



SO 



^so 



be, — d 1 



(21) 

where dl — \d\ for counterclockwise rotation of the elec- 
tron and dl — —\d\ otherwise, dr is the projection of a 
on the radius vector. Eabc,— d is the oriented area of the 
trapezium formed by the vectors a, 6, c and —d. It is 
positive if the contour abc, — d is positively oriented and 
negative otherwise. The last term in (21) takes into ac- 
count the finite width of the ring. It leads to a phase pro- 
portional to the area embraced by the trajectory, anal- 
ogously to the finite width effect on the Aharonov-Bom 
phase^**. 
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Formulae (19) and (21) yield the following dynamical 
equation for the spin evolution operator 



Here 
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■■ a/Lso , 



hit), 
a 



(22) 



(23) 



(24) 



a is the ring radius (distance from the center to the mid- 
dle line of the ring), and cp is the time derivative of the 
polar angle (p counted from the negative direction of the 
OX-axis. The deviation h{t) of the electron trajectory 
from the middle line (see Fig. 2) is taken positive for 
points outside the circle formed by this line and negative 
otherwise. 

Dynamical equation (22) simplifies if we change the 
system of coordinates. First we perform transformation 

to the system OXYZ rotating together with the electron 
(OZ-axis coincides with the OZ-axis, OX and OY axes 
rotate with the angular velocity (p around OZ-axis). In 
this system the evolution operator takes the form 



S = R,{t)SR-\0) , 
R,it) = exp{ia,(l){t)/2} . 



(25) 
(26) 



Using (22) and (25) and taking into account that 
R~^{0) = 1 one can verify that S satisfies the equation 

^ = |-^^-^r/a.|5, (27) 

where ae is a projection of & onto the vector 

= {27/C,O,l/C} , (28) 

C=\/r+4^, (29) 

and 



The evolution of S' is determined by 



dS' _ .&n~ 

dt ~^ 2 



-Wo 



277 



(33) 
(34) 



The physical meaning of the TR system is rather trans- 
parent. At vanishing ring width and zero magnetic field 
this is the coordinate system where the effective mag- 
netic field produced by the Rashba SOI is parallel to the 
OZ'-axis. 



Evolution of spin state along a trajectory in 
terms of polarization vector 



Transition from 5 in (11) to 5' yields 



(35) 



where 



RyRzO'iR'^Ry 



(36) 

RyR,(Tj{t)RtR+ , (37) 



Ry(fjRy ^ 



SajS+ . 



It is easy to sec that fj^' in (36) is an image of ai af- 
ter transformation to the TR coordinates. Analogously, 

[(7j(t)] is is an image of (Jj{t) after the same transfor- 
mation. Note, that both 5- and [aj (t)] depend on time: 
al — due to rotation of the TR system of coordinates, 

[(Tj(f)] due to rotation of the electron spin (described 
by the evolution operator) and rotation of the TR sys- 
tem of coordinates. Further, both ct- and [crj (i)] can be 
decomposed into the sums over ap, 



1^0 = -C<P ■ 



(30) 



Next, we rotate OXYZ around OF-axis,until OZ 
comes parallel to 0. The new coordinate system denoted 

as OX'YZ' will be called below the tilted rotating (TR) 
system. In TR coordinates the evolution operator takes 
the form 



— S R'y — ^ ^^y 

Ry = exp{iayXa/2} . 



(31) 



The angle xo between OZ and OZ' is defined by the 
relations 



cosxo = 1/C, sinxo = 27/C. 



p=i 
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p=i 



(38) 
(39) 



The superscript "0" at Ap, in (38) indicates, that Api 
transforms into A^- with ago — > 0. The unit matrix ctq 
is not present in these sums because the traces of and 

[(Tj(t)] arc zero, as can be seen from (36) and (37). Using 
(B2) from Appendix B one can check, that Apj are TR 
coordinates of the electron spin, provided that the initial 
spin was Sj. Analogously, A^^ arc coordinates of in 
(32) the TR system. After substitution into (35) the sums 
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(38) and (39) lead to Kij written in the form of scalar 
product, 



It is convenient to represent the angle $ in the form $ = 

ip + d, where 



aN{a;„a°„a§j 

A, ={Ay,A2,,A3,} 



(40) 



(44) 



The components of A? in (40) being the coordinates of 
ei in the TR system are defined only by its orientation 
with respect to the original system, that is by the angles 
(j) and xo- The angle cj) at the end of the trajectory is 
(f) = 7r + 27rn, where n = ±1, ±2... is the winding number. 
The angle xo is also fixed, sec (32). Hence, components 
of A? at the end of a trajectory are constants defined 
explicitly by Eq. (36): 



A" 



-io-^ 
C ' C 

AO = {0,-1,0} , 



A" = 



27 1 



(41a) 
(41b) 
(41c) 



With constant A° averaging of Kij over trajectories re- 
duces to averaging of A^-. We shall perform this aver- 
aging with the use of the distribution function V of Aj . 
This function will be obtained from the Fokker-Planck 
equation which will be derived and solved below. 

The equation of motion for [crj (t)] is found from (37) 
and (33). In view of the expansion (39) it can be written 
as an equation of motion for A^ : 



dAj 

dt 



(42) 



This equation should be supplemented with the initial 
conditions 

A.|,=o = {l/C,0, 27/C} , (43a) 
A,U = {0,1,0}, (43b) 
A.U = {-27/C,0, 1/C} , (43c) 

which can be derived from (37) and (39). 



For these new variables we obtain the system of equa- 
tions: 



e =27(^i^+WM/(t))sin(V' + (5) 
S = uw{t)^ [1 + 27C0s(^ + S) cot e] 



(45) 



where we have introduced the frequency u!w associated 
with the finite width of the ring, 

cjwit) = f^it)uo{t), Kt) = (^^^ ' ^ . (46) 

At the weak enough magnetic field ujh "C ojq- If, in 
addition, the ring is narrow, h{t) <C a, then /Lt(t) <C 1 
and LJw ^ ^0, as follows from (46). Returning to (45) 
we thus see that & and 6 are "slow" variables, while ip isa, 
"fast" variable. Such a separation of variables simplifies 
considerably the Fokker-Planck equation, which will be 
derived below from the system of stochastic difi'erential 
equations (45). 



D. Derivation of the Fokker-Planck equation 

Parameters of the Fokker-Planck equation for the 
distribution function V{Q,6,ip) are determined by the 
drift Aq, As, A^ and diffusion Bij;i,j = {Q,5,ip} 
coefficients^^. Let us first consider the difl[usion coeffi- 
cient B^^ given by 



B^^ = lim 



(AV)^ 



At^O At 



(47) 



where Aip is the increment of the stochastic process tp{t). 
It follows from Eq. (44) that {Atpy = CJA^. On 
the other hand, it was shown in Ref. 15 that the wind- 
ing number w of the classically chaotic trajectories has a 
Gaussian distribution 



Stochastic differential equations for the angles 
determining the position of an electron and 
direction of its spin polarization 



It follows from (39) that A^ = Tr <^ [aj{t)] W2 = 1. 



This allows to describe Aj by only two variables, the 
polar angle $ and azimuthal angle 6. Eq.(42) is then 
reduced to 



e 



^77 sin $ 



$ = wo(i) + ^r][i + 27 cos $ cot 9] . 



V{w\T) = (27rT/Ti) 



-w 



-1/2 / 

^"Pi2T77T 



(48) 



where the constant Ti is the characteristic time of one 
turn, {w^{Ti)) = 1. This equation means that the wind- 
ing of trajectories is a diffusion process. One can extend 
(48) to a range of w = (j)/2TT < 1 assuming that a par- 
ticle advances diffusively along a ring arrn. Such situ- 
ation takes place if the rotation direction (l)/\(t>\ changes 
many times, while passing the angular distance < 27r. 
Moreover, the angular distance between two consecutive 
changes of a rotation direction must be small enough to 
ensure a small change of $ and tp. Hence, as follows from 
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Eq. (44) this distance must be <C C~^- Wc assume that 
scattering from ring boundaries and spatially fluctuating 
potential make this condition satisfied. 

From (47) and (48) with w = ^/27r we immediately 
find 



B. 



lim 



(AV)^ (27rC)^ 



At-yO At 



Ti 



(49) 



where the overline denotes averaging over trajectories. 
Other diffusion coefficients, as well as drift coefffcients, 
are conveniently expressed in terms of the diffusion coef- 
ficient of the auxiliary stochastic processes u{t), which is 
defined by its stochastic differential 



du = LOw{t)dt . 



(50) 



An assumption that the deviation h{t) of a trajectory 

from the middle line of the ring and the angle (/> fluc- 
tuate independently leads to the absence of correlations 
between the stochastic processes ^p{t) and u{t). Assum- 
ing also a uniform distribution of h over the width d of 
the ring arms, we find from Eq. (46) 



lim 

At-)-0 



A'^ Am 



At 



0, 



Buu = lim 

At->-0 



{AuT 
At 



12a? 



B. 



(51) 
(52) 



Further, Eq. (45) can be used to express the small 
increments AQ, Ad, Atj; in the form of integrals over 
time interval At. Then, after averaging procedure the 
limits At ^ must be taken. For example, Aq = 
limAt^o Ae/At, Bqs = limAt^-o AQA5/At. We thus 
obtain 



Ae 
A^ 
As 



= uoH^ sin {,p + 5) + ^2-i cos (V' + (5) F , 
= 0, 

= - ^27 sin (V- + 6) [cotQ 

+ 27^1^cos(V' + (5) 



= sin2(^ + (5)(27)2i?„ 

-'-'7/1/ 1 



Bee 

Bss = 

Bes = 2j sin {ij + S)FB^ 
B^Q = B^s = , 



(53) 



(54) 



where 



F =1 + 27cos(V' + 6) cot e . 



These coefficients should be inserted in the Fokker- 
Planck equation 



d_ 

di 



{VsmQ) = 



where the variables Xi [i = 1,2,3) denote 0, 5, tp. The 
probability density V is normalized in such a way that 
an integral of PsinS over 9,(5 and ij) is 1. In this way 
we arrive at the following equation 



sin 9 



&P B. 



dt 



1 l + cos[2{tP + S)], 



2-i f dV 

-ujH-^sin{il) + 5) ( — sin9-FPcos9 \ + .. 



(55) 



where the dots stand for other terms that are propor- 
tional to Buu or ujh- All the terms in r.h.s, except for 
the last one, do not contain derivatives over ijj. For a 
narrow enough ring and weak magnetic field, and 
ujh ^ B^,^. So, diffusion in the space of the "slow" vari- 
ables 9 and 5 is indeed slower than diffusion through the 
"fast" variable V- Hence, the diffusion equation can be 
averaged over the fast variable. After averaging (55) over 
xl) we arrive at the Fokker-Planck equation of the form 



dt C 95 ^ ""^ sin 9 99 

-+7W9U^ 



Bi, 



""®a9j 

' 2 dip"^ 



(56) 



This equation should be solved together with the initial 

conditions for each of the vectors Aj . 

r{Q, tjj, S\t = 0) = (5(cos 9 - cos el)S{S - Sl)S{tjj) , (57) 

where the angles 9q, 6q defining the initial positions of 
the vectors Aj are found from (43). Using Eq.(32) ^0 
and So can be expressed in terms of xo- The initial value 
of ip should be zero since it is proportional to the initial 
value of 6, which is zero. 



E. Solution of the Focker-Plank equation 

After Laplace transformation with respect to time and 

Fourier transformation with respect to S and ip the equa- 
tion (56) is reduced to the ordinary differential equation 



Xv 



where 



1 d ^ d \ 

sm9— — ^ — > V 



sinQdey dej sin^e, 

= (5 (cos 9 -cos 9^) , (58) 



v = j^B^ue'^'iV., 



V = V{&, K, q\p) = j d^dS e-^(«''+9'^) 
r+00 

X / dte-P*V{Q,S,ilj\t) 

_P + IK^ + K^Buu - 7') + 
~ -72 R 
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The solution of Eq. (58) can be expressed in terms of 
eigenfunctions of the linear operator in l.h.s. of this equa- 
tion. In our case they are the associated Legendre func- 
tions p1'*'(cos9) and we obtain 



E 



2n 



+ 00 



X E \^M!:l:li; ^-'(cos9)ii-'(cose^„) 



2 {n+\K\)\ 

X e-B„„[l«r(i-7^)+7Vn+l)]*^ (59) 



In this equation only a factor in front of the first 
sum depends on ip. It is clear that this factor deter- 
mines a probability distribution of V'- At the end of tra- 
jectories (at the exit from the ring), when (f) = 2ttw, 
w = ±1/2, ±3/2..., it coincides with the winding number 
distribution (48) . The remaining part of (59) is evidently 
the conditional (for given ip) probability distribution of 
S and O. This function can be used for averaging of 
the polarization vectors Aj over trajectories with a given 
winding number w. 



F. Averaging of the polarization vectors Aj 

Since the unit vectors Aj in Eq. (42) are defined by 
their respective polar and azimuthal angles 6 and one 
can calculate easily their average values using Eq. (59) 
and taking into account that ^ = tp + S. We thus arrive 
at the following expressions for the averages (Aj)!^^^^ at 
fixed trajectory duration T and a given ip (winding num- 
ber), 

(Ax,)|v.,T = sine^o cos (^V + + S^j e"^/^^ , (60a) 

{Ayj)\^,T = sine^o sin (^^ ± ± d^j e"^/"^ , (60b) 

(60c) 



(A.j)|^,T=cose^oe-^/T. 



The parameters t± and T|| have been introduced to char- 
acterize relaxation rates of the electron polarization due 
to the finite width of the ring, 



1 _ 1 + 



TJ_ 4 

- = 27'S, 



(61) 



As follows from Eq. (60), perpendicular to the OZ'-axis 
components of Aj decay with the rate l/rj_, while the 
decay rate of parallel components is 1/t||. Wc recall that 

in the rotating system the direction of the OZ'-axis is de- 
termined by the vector 6, see (28). At t = the rotating 

system coincides with the original one. Hence, the elec- 
tron polarization in the ring relaxes with the rate t± (rii ) , 



if the polarization of the left reservoir is perpendicular 
(parallel) to 6. Besides relaxation associated with finite- 
ness of the ring width, there is an additional relaxation 
channel due to the magnetic field, that will be discussed 
below. 

To complete calculation of the spin conductance given 
by Eqs. (14) and (40), the scalar products A? • (Aj)|^^T, 
where A? are given by Eq. (41), must be aver- 
aged over tjj and T. Averaging over ip is performed 
with the use of (48), by substituting w = —ip/27r(. 
As for the distribution over T, in the case of classi- 
cally chaotic systems one should use the exponential 
functioni6 ■p(T)=T-i exp {-(T - Tq) /t}, where r is the 
mean escape time of a particle and Tq is the shortest 
trajectory duration. The results of calculation for gyy, 
as well as for polarization rotation angles in a magnetic 
field are shown in Fig. 3 and Fig. 6 respectively. It is 
seen from these plots that the AC oscillations magnitude 
and spin rotation angle strongly depend on the parame- 
ter Ti/t + Ti /t± , which controls the trajectory winding 
number during the particle spin lifetime. If this parame- 
ter is small, w is large and AC oscillations are strong. In 
this regime one can write simple analytic expressions for 
tensor components of the spin conductance: 



(5xx) =-50^{Q + 472m} 

{9yy) = -90 Q , 
(5zz)=5o^{47'Q + M} , 

idxy) = -{gyx) = 9o\r, 



C 

27 



(62a) 
(62b) 
(62c) 

(62d) 



{9.z) = -{gz.)=go^{Q-M} , (62e) 



27 

{9yz) = {9zy) = 90-^R, 

where Q, R and M are given by 



Q = costtC 

R = costtC 
1 

M = 



1 + ^ + l^sinVC 



(l + ^ + l^sin^^C) +(^) 



1 + 



(62f) 

2 ' (63a) 
2, (63b) 
(63c) 



Note, that expressions (62) were obtained under an as- 
sumption that the magnetic field is not too strong, so 
that WffTi/C < y/6Ti{l/T+l/T±), while Ti/2To > 1. 



V. RESULTS AND DISCUSSION 

Wc start our discussion from the analysis of the finite 
width effects. For simplicity, we will consider the {gyy) 
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g 



yy 
go 




FIG. 3: Attenuation of spin current oscillations due to finite 
widths of tlie ring and leads, at a zero magnetic field, gyy 
is the j/-component of the spin conductance along the y-nxis, 
see (3) in the text, go is the transmitted flux of electrons 
per unit energy interval, per one spin projection of the inci- 
dent electrons. The variable is expressed in terms of the 
spin-orbit length Lso sls — ^\ -f 4{a/ Lso)'^ ■ Three curves 
correspond to parameter sets 

(a) a = 1 /I , d — 150 nm , W — 150 nm , g = 1/2 ; 

(b) a = 1, 5 A* , d = 300 nm , W ^ 100 nm , g = 1/5 ; 

(c) a = 1, 5 /i , d = 300 nm , W = 30 nm , q = 1/5 . 



component of the spin conductance matrix in the regime 
of large winding numbers when analytic expressions (62) 
are valid. In the absence of the magnetic field the spin 
conductance is given by 



idyy) = -5oCOS7rC 



^ sin^ ttC + i + 



(64) 



The denominator in (64) gives rise to a set of peaks 
with maxima at ^ = Cm = m, where m is integer. It is 
easily seen that for T^^ ^ -I- rj^ the peak's broad- 
ening Ag <C 1. Due to this inequality Eq. (64) can be 
written in the vicinity of peaks in a more simple form: 



90- 



Ti 
2tt^ 



[ACf 



Ti 



1 



(65) 



where AC — C, — C,m- From Eq. (65) is expressed as 



/5 













(66) 



Note, that as follows from (C6), t±_ depends on C- Hence, 
Ag depends on the resonance Cm position. 

In the case of a long particle lifetime r ^ one ob- 




B 



FIG. 4: A schematic picture explaining the effect of finite 
width. Since the element of the trajectory passed forward 
and backward does not give a contribution to the evolution 
operator, one can replace the trajectory 2 by 2'. After that 
it is obvious that the difference between 1 and 2 is that the 
latter contains the loop ABCD passed counterclockwise. 



tains from Eq. (C6) 



d 8(1 + C^) 



3Cn 



(67) 



For example, the broadening of the third peak (Cm = 
3,7m = a/S) is Ag ~ iAd/a. It is a quite noticeable 
value for a typical ratio d/a ^0.1 

So, the first obvious effect of the finite width is the 
broadening of the spin current oscillation peaks. The 
physical origin of this effect is the increased relaxation 
rate of the spin polarization. This relaxation is caused 
by incoherent superposition of polarizations coming from 
the trajectories encircling slightly different areas in a ring 
of finite width. The situation is elucidated in Fig. 4. 
This picture shows that the finite width results in adding 
random loops breaking the coherency of the trajectories. 

In addition to the broadening, the increased relaxation 
rate leads, evidently, to a reduction of the peak intensity. 
This is explicitly given by 



(5: 



yy I Q=Cr, 



go 



1 



(68) 



which is the spin current magnitude exactly at maxima 
(minima). From Appendix C the ratio of times in the 
denominator of (68) can be expressed as 



r ^ 1 7^(1 + C') a 
r_L 6nq C^ ^ 



2a 



(69) 



This expression shows that the finite width effect is sup- 
pressed fast with smaller d/a. 

Now let us focus on magnetic field effects. The first 
effect is that the components (gyz), (gzy), (gxy), (gyx) of 
the spin conductance are no longer zero. One can verify 
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FIG. 5: Narrow loops of the trajectory which break the co- 
herence of the trajectories in the presence of the magnetic 
field. 



from Eqs. (14), (11) and (8) that they appear because 
the reflection symmetry with respect to the XOZ-p\ane 
is broken by the magnetic field. The physical meaning of 
such nondiagonal components can be explained in terms 
of the effective polarization Peff{ej) on the exit from the 
ring, see Eq. (16). For example, nonzero (gzy) and (gxy) 
are associated with a rotation of Peff(ey) with respect to 
the polarization By of the left reservoir. It is convenient 
to consider a projection of Peff{ey) onto FO^'-plane. 
Then, the rotation angle 0p of this projection can be 
calculated from Eqs. (62) and (63). For the m-th peak 
this angle is given by 



tan( 



27 



Cm 



(70) 



The nonzero {gyz), idyx) components can be interpreted 
in a similar way. We note, that due to the linear depen- 
dence on luh, the sign of 0p changes together with the 
magnetic field. 

Another effect of the magnetic field is a reduction of 
the spin current. Let us consider a trajectory which con- 
tains a narrow loop, see Fig. 5. If the magnetic field is 
ignored and only the SOI effect is taken into account, af- 
ter passing the loop the polarization P does not change. 
That is because on the upper and the lower parts of the 
loop P rotates in opposite directions, according to op- 
posite directions and of the SOI fields, see Eq. 
(8). A magnetic field, however, causes rotations of P in 
the same directions. Hence, evolutions of P along tra- 
jectories with and without the loop becomes different. 
This introduces an additional decoherence leading to the 
spin current reduction and broadening of its oscillation 
peaks. Using (62) and (63) one can derive the following 
expression for the magnitude of the effective polarization 




FIG. 6: Rotation of the effective polarization vector Peff{e.y) 
on the exit from the ring in the presence of the magnetic 
field, (jip is the angle between the initial polarization By of 
the electrons in the left reservoir and projection of Peff{e.y) 
onto the YOZ plane. The magnetic field strength is lOOG. 
Curves (a), (b) and (c) correspond to the same parameter sets 
as in Fig. 3. 



exactly at maxima (minima) 

= %^\l\{9yy)c=U? + \{9zy)c=cS + \{9xy)c=cJ'' 



1 + 



-1-1/2 



(71) 



provided that the ring is narrow enough, t/tx <C 

Let us consider a dependence of the spin conductance 
on C in the presence of the magnetic field. In the the 
vicinity of the m-th peak, instead of (65) wehave 



{9yy) = 5o(-l) 



1 



^(AC)2 



i + _L + ^(AC)2] +(^) 



(72) 

As can be seen from this equation, the relaxation mech- 
anism associated with the magnetic field gives rise to an 
additional broadening of spin current peaks. Their width 
can be evaluated from Eq. (72) as 



2Ti 



1 1 

T rj_ 



UJH_ 

Cm 



1/4 



(73) 



The discussed above effects are determined by char- 
acteristic times r, Ti and tj^. These times have been 
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evaluated in Appendix C using a simple model of scat- 
tering from a bumpy ring boundary. Fig. 3 and Fig. 6 
demonstrate the effects of the ring width and magnetic 
field on behavior of the spin conductance as a function 
of C = ^/l + ^a/Lso)^. We took kp =i 2.5 x 10*^ cm"! 
and Lso varying in a wide range. In IriAs based quantum 
wells SOI can be quite strong with Lg,, being as small, 
as ~ 100 nm^^'^'^. For o = 1/x this gives ( = 20. Curve 
(a) in these figures corresponds to the escape time much 
shorter than Ti and t± . The winding number is not large 
and AC resonances are broad. There are no noticeable 
effects associated with the finite ring width. Also, the 
magnetic field effect is relatively weak. The width effects 
are seen on the curve (b), in Fig. 3. A reduction of the 
oscillation amplitude seen in the figure is in a qualitative 
agreement with Eq. (68), although for considered param- 
eters this equation can not be fully applied, because the 
winding number is not large enough. The winding num- 
ber is larger for the third set of parameters, (c). The 
finite width effect becomes stronger , leading to a faster 
decreasing of the oscillation amplitude. Also stronger is 
the magnetic field effect on polarization rotation in the 
a;j/-plane (see Fig. 6). Strictly speaking, our semiclassi- 
cal theory can not be applied to this case, because the 
number of propagating channels in leads is not large. We, 
nevertheless show this result in order to demonstrate a 
trend: for reasonable ring sizes the regime of large wind- 
ings with sharp AC resonances can be achieved only at 
small lead widths, or by means of barriers between leads 
and the ring, resulting in the long r. 

We note that the magnetic field effect on polarization 
rotation is rather noticeable even at relatively weak 100 
Gauss magnetic fields, as can be seen at Fig. 6. For 
example, for C = 1.6 and parameters (c), the rotation 
angle can be as large as 24°. 
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Appendix A: Landauer formula for the spin current 

Since the spin current density is an additive one- 
particle dynamical observable, its average value at the 
point R and t = is given by 

( Jy (R, t = 0)) = Tr { A (i = 0)J[^ (R) } (Al) 

Here the one-particle distribution function /i describes 
the open system consisting of the leads and the ring. The 
one-particle operator 

j,^^^^virjm + Mm (A2) 



represents l-th component of the current density with 
spins polarized along j-th coordinate axis, vi is ^-th com- 
ponent of the electron velocity operator. Since we calcu- 
late the spin current in the asymptotic region of the right 
lead, where the magnetic field and SOI are zero, the op- 
erator -Oj may be written simply as pi/m*, where p is the 
electron momentum operator. The polarization density 
Vj(R.) is defined by 

Vj(R) = pn<fj, (A3) 

where /5r is the density operator. In the the coordinate 
representation the latter is given by 

PR = 6{r - R) . (A4) 

It should be noted that Eq. (Al) represents a polariza- 
tion current density, rather then the spin current density, 
which is twice smaller. For convenience we will use, how- 
ever, the latter name. 

For noninteracting electrons the evolution of the dis- 
tribution function is described by the equation 

iH^ = [H,h] (A5) 

where H is the one-particle Hamiltonian for the system 
"leads-hring" . A formal solution of Eq.(A5) may be writ- 
ten in the form 

hit) = Uit,to)fi{to)U+it,to), (A6) 
U{t,to) = e-'"^*-'°^/''. (A7) 

Further, we take into account that the system \mdcr con- 
sideration has an asymptotic region where an electron 
is effectively decoupled from the ring. In this case the 
methods of the scattering theory may be applied directly, 
without adiabatic switching off the scattering potential 
a,t t = ±oo. First, let us write down /i(0) in the form 

/i(0) = lim f7(0, -r)f/o+(o, -r)[/o(0, -T) 

T— >-+oo 

X h{-T)U+{0, -T)f/o(0, -T)f/+(0, -T) , (A8) 

where the unperturbed evolution operator Uq is obtained 
from ( A7) by replacing H with the " unperturbed" Hamil- 
tonian Hq. The latter is obtained by removing the ring 
and elongating the leads to meet each other. In Eq. (AS) 
one easily recognizes the familiar MoUer operator r2+ of 
the scattering theory^*, 

n+= lim U{0,-T)U+{0,-T). (A9) 

T^-l-oo 

This operator maps the wave function {ipin) describing a 
particle state at i = in the absence of the ring onto the 
actual state \tp{t = 0)) : 

\^{t = 0)) = n+ ■ (Alo) 
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From (A9) and (A8) we obtain 

/i(0)=O+/r(f2+)+ , (All) 
where, by analogy with \tpin), the function 

/r= Um Uo{0,-T)h{-T)U+{0,-T) (A12) 

I — >+oo 

can be interpreted as a distribution function of the sys- 
tem at t = in the absence of the ring. The trace in 
(Al) can now be rewritten as 



(jy(R,t = o)) = ^(*ii/rK2) 



X {i2\ (0+)+ (R) n+ In) . (A13) 

Since the unperturbed problem does not involve SO in- 
teraction, a convenient choice of the basis vectors \i) in 
(A13) is 



\i) = \b} ® |a) , 



(AM) 



where the eigenvector \b) corresponding to the unper- 
turbed Hamiltonian describes the electron orbital motion 
and I a) is the eigenvector of az corresponding to its eigen- 
value a. Further, the slab geometry of the unperturbed 
problem suggests that \b) is taken in the form 



\b) = |fc) (g) \p) (g) |m) 



(A15) 



with the eigenvectors |A;),|p),|m) describing a particle 
motion along OX, OY, OZ axes in the absence of the ring. 
Hence, the corresponding wave functions are 



Wk{x) = {x\ k) = ■y^e*'=^ , 

'"p{y) = {y\p) = \It- ^Mkyv) 



ky = —p, {p= 1,2,...) 



(A16) 

(A17) 
(A18) 



and similarly for the wave function Um{z), m = 1.2,..., 
in z-direction. We took periodic boundary conditions in 
x-direction, where L is the total length of the system. At 
the slab interfaces the wave functions Vp{y) and Um,{z) 
satisfy the hard wall boundary conditions. 

Unit vectors parallel to polarizations of the left and 
right reservoirs will be denoted as and i/^, respec- 
tively. Accordingly, we define the operators g^l.r = 
J2i^i'^i'^ with eigenvectors \i/^'^a) corresponding to 
polarization projections a =t,4 onto i/^ and . Since 
particles with different spins are distributed in reser- 
voirs according to their respective Fermi distributions, 
the magnitudes of the reservoirs polarizations arc deter- 
mined by the differences 5ijll,r = Ati/^.«t~^i^^'"4- of chem- 
ical potentials of spin up and spin down (relative to v^'^) 
electron gas components. Therefore, assuming that the 
unperturbed distributions of particles moving to the right 



{k > 0) and to the left (fc < 0) are given by the Fermi 
distributions in the left and right reservoirs, respectively, 
we can write 

/r = Y.®{k) n^,.., (|6) (6|) ® {\v^a) {v^a\) 

+ ^ e (-fc) n(,,.H, (|6) (6|) (|i/«c7) {y^a\) , (A19) 

h,cr 

where ©(•) is the Heaviside step function and 
'nb,v^a = ifJ-i^^a - E) , ni,^„R„ = e {n„Ra - E) (A20) 

are the Fermi distributions in the left and right reservoirs 
for particles with the energy E. Note, that Eq. (A19) 
was written under the assumption that contacts between 
reservoirs and leads are adiabatic (no scattering from the 
contacts). 

We assume for the average chemical potentials hl = 
fjLR, where Hl/r = {ij,^l/r^ + fii,L/R^_)/2. So, the chemical 
potentials of unpolarizcd reservoirs coincide. Denoting 
them we write the distribution function correspond- 
ing to the unpolarized reservoirs in the form 



(A21) 



where = Q{pu — E). Evidently, does not give any 
contribution to {Jij{'R,t = 0)). Therefore, it is conve- 
nient to subtract this function from /™: 

/r = /T + + ,5/« (A22) 

^fi =J2®ik)6n,^,., (|6) (6|) ® (l^^^a) (i.a^|) (A23) 



^fi =Y.e{-k)Sn,^,R, (|6) {b\) ® {\u^a) (u^a\) 

b,a 

(A24) 

Sub^^L.R^ = Ub^^L.R^ - (A25) 



Denoting corresponding contributions of 6/^ and Sf^^ to 
the spin current as ( Jy (R, t = 0))l and ( Jy (R, t = 0))/?, 

we arrive at 

(Jy(R,t = 0)) = {Jij{R,t = 0))L + {Jij{'R,t = 0))R. 

(A26) 

The projectors ji^-^cr) (z/-^(t| and \i'^(T)(^i'-^'a\ in (A23) 
and (A24) can be expressed in terms of the Pauli ma- 
trices a^L,R and the unit matrix ao using easily verified 
relations 

Straightforward calculations then give 



{Jlj{R,t = 0))L,R = 



E 



e{±k)-^ 



" ..L.R 



ML' 2 — <-E</Jr;H 2^ 



X (6a2| {VL+Y (R) n+ \bai) , (A27) 
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where \ba1_2) = |&) ® and upper (lower) sign in 

the argument of 0-function corresponds to the index "L" 
("R"). The leads are assumed to be thin enough in z- 
direction, so that only the levels \b) = \kpm) with m = 1 
are occupied and contribute to the sum in Eq. (A27). 
For simplicity, we denote 



where 



\kp) = \kp, m = 1) 



(A28) 



It is convenient to change in (A27) the summation over 
k by integration over E. To do this, we introduce the 
vectors 



\Ep)^^^ = ^f^)\kp) 



(A29) 



where ± signs relate to A; > and A; < 0. t'p(-B) is the 
one-dimen sional density of states in the p-th channel, and 
Vp{E) = ^2{E- Ep)/m* is the electron velocity in p-th 
channel characterized by the kinetic energy in y-direction 
Ep. Using definitions (A29) one can write 

E ^i±k){b\.\b} 

b 

IJ.U j^<£;<jUf7H 



= E 



dE ^^\Ep\ ■ \Ep)^^^ , (A30) 



Further, in the limit L +00 considering E as a. contin- 
uous variable one gets the normalization condition 



^^\E'p'\ Ep)^^^ = 5p.p6{E' - E) . 



(A31) 



Using this condition and substituting (A30) into 
Eq. (A27) we find 



(J,,(R,i = 0)KH= E / ^ '^^^ 

X W(Epa2| J{' (R) n+ \Epai)^^^ , (A32) 

{\Epai,2)^^^ = \Ep)^^^ (g) |ai,2)). In its turn, the total 
spin current through the cross section of the right lead is 
given by 

{J^{X,t = 0))L,R = J dYdZ{J^,{K,t = 0))L,R. 

Using Eq. (A32) we obtain 

{Jj{X,t = 0))L,R 

= y / dE^ dYdZ 

p,ai,a2 2 — •' 

x^^\Epa2,+\J^' {K)\Epau+)'^^^ , (A33) 



\Epa„+)^^^ =ft+\Epa,)^^\ s = l,2. 



The vectors \Ep as, +)^^^ are known^^ as the scatter- 
ing states associated with the "in" asymptotes ("inci- 
dent waves") \Epas)^^K Since the point R in the r.h.s 
of (A33) is located in the asymptotic region of the right 
lead, only the asymptotic behavior of the wave functions 
= SP"2,+)'^^ and <^g,\^^^^+(R/3) af- 
fects the calculation of the matrix elements in (A33). 
Thus, we may write ?5'^^q,2,+ (RQ!) and (j^'^E^ai^+i'^'^) 
the sum of transmitted waves while <l)'Epa2 +(I^Q;) and 
'^Bp^c«i,+ (^'^) ^ ^^'^ °^ incident and reflected waves. 



'„^_+(Ra) = E 4'p'{m^Ei'a''i^<^) (A34a) 



+ E 4'TiE)^Ei"c."{^(^l (A34b) 

p" a" 

where t"y,p^{E) and r'^„p^{E) denote transmission and 
reflection amplitudes, respectively. According to (A29) 
and (A28) 

C''a"(R«) = ^l'^Ux)up„{Y)v,{Z)xc"{a) , (A35) 
where 



Vpn{E)w±k{X) = 



^^-Khopii {E) 

(A36) 

and Xa"{ct) = ct")- k in (A36) is the positive solution 

of the equation E — Ep" = fi^k"^ /2m* . Note, that due 
to our choice of the prefactor in (A29), the functions 
^p"s(X) "carry" the same flux (27r?i)~^, independent on 
the channel number p" . As a result, the transmission and 
reflection amplitudes satisfy the flux conservation law 



p a 





2 






1 La"a2 


+ 






1 Vp 


P"P 





Using Eq. (A2) and Eqs. (A34) - (A35), we transform 
Eq. (A33) into 



{Jj{X,t = 0))L 
_ 5jlL_ \- 1 



p,p' 

{Jj{X,t = 0))R 



2-Kh 



NmilJ'u)-^ Tr{ijj,BCTj} 



(A37b) 
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where tpip{fiu) is the 2x2 matrix composed of the 



transmission amphtudcs t'^Jp{fMj),'tlfp{fMj) , 



tpfpdiu) ■ The matrix fp'p{iiu) is composed in a simi- 



lar way. Taking into account that a^, 
denoting 



a ■ V = GiVi and 



9% = 



4 



2^ E {t+>(M£/) ip,p{iiu) ^i} (A38a) 
^Tr \^f^,p{^iu) ^ j rp'pinu) ^ i] 



1 

2^ 



-Nm{m)k3 

(A38b) 



we rewrite (A37) in the form 



(J,(X,^ = 0))i,i^ = V 5 



(A39) 



We shall call the 3x3 matrices and as spin conduc- 
tances. They determine the response of the spin current 
to the polarization of the left and right reservoirs, re- 
spectively. For example, if the left reservoir is polarized 
along the s-th axis, i^^ = Bg (Cg is the s-th coordinate 
ort), and the right reservoir is unpolarized, then accord- 
ing to (A26) and (A39) {Jj {X,t = 0)) = .gjlj^Mi and gf^ 
proves be the proportionality coefficient between Sfx^ and 
j-th component of the spin current. 

It is convenient to express the quantities v^Siil and 
vf^S^iR in Eq. (A39) in terms of 2D spin polarization 
densities {Vi{XY)) ^ and {Vi{XY)) n corresponding to 
chemical potentials of the left and right reservoirs. At 
small SfiL and djiR <C Ep we have 

{V^{XY))L/R - Npyl'^'SiiL/R , (A40) 

where Np = 27rm*//i^ is the 2D electron state density. 
Substituting this expression into Eq. (A39) we rewrite 
the latter in the form 



(J, (X, t = Q))l,r = ^ E4'''(^^(^^)) 



L,R ■ 



(A41) 



Combining (A26), (A39), and (A41) we obtain finally 

(J, (X, i = 0)) = (J, (X, t = 0))l + (J, (X, t = 0))k 



Yg^,{r.{XY))L+Ygf,{r.{XY))R \ , 

1 i i ) 

(A42) 



This expression gives the spin current in the right lead 
in terms of the left and right reservoir polarizations. 

Comparing our expression (A38) for the spin conduc- 
tance with the Landauer formula we see that they are 



different in the way that our spin conductance is written 
in terms of both transmission and reflection coefficients. 
This difference is of principal character since due to non- 
conservation of spin current we can not express the con- 
tribution {Jj {X, t = Q))r (containing the reflection coef- 
ficients) via the spin current in the left lead transmitted 
from the right. 



Appendix B: Evolution of polarization along s-th 
trajectory 

If we define the polarization of an electron as 

P = Tr{pd-}, (Bl) 

then, according to Ref. 29 the density matrix describing 
its spin state may be represented in the form 



(B2) 



where tfo is the 2x2 unit matrix. After passing s-th 
trajectory, the spin density matrix is transformed into 



SspSf 



(B3) 



where Ss is the spin evolution operator given by Eq. (8). 
Substituting (B2) into (B3) and then (B3) into (Bl) we 
obtain the polarization of the electron at the end of the 
s-th trajectory 



Writing down this vector equation in components we ob- 
tain 

pi — 1{S p. 

where K-^ is given by Eq. (11). Therefore, we see that 
Kfj gives i-th component of the electron polarization at 
the end of s-th trajectory, provided that at its beginning 
the electron had the unit polarization along the j-th axis. 



Appendix C: Characteristic times 

In this Appendix we derive expressions for the particle 
lifetime r, the relaxation time associated with the finite 
width of the ring t± and the characteristic time of one 
turn Ti. The widths of the ring (d) and leads {W) will 
be assumed to be much less than the radius a. 

Let us start with r. This time is determined by the 
shortest of two times: the mean time of particle escape 
from the ring and the dcphasing time associated with in- 
elastic electron-electron and electron-phonon collisions. 
We will assume that the temperature is low enough to 
neglect the latter effect and will focus on the escape time. 
Any electron trajectory inside the ring is a set of straight 
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segments. The probability that a current trajectory seg- 
ment is the last one before escaping from the ring is r^/r, 
where t/j is a time interval between two consecutive col- 
lisions with ring boundaries. On the other hand, the 
same probability may be written as 2W/A-Ka. Thence, 
r ~ r/i(27ra/W). In its turn, th can be estimated as 



d_ 

Vf 



1 



cos^ 



(CI) 



where 6 is the angle between the particle velocity and the 
radius- vector. Eq. (CI) is valid for 6 not too close to 7r/2, 
namely \tt/2 — ^| ^ yj2d/ a. The average in this equa- 
tion is calculated assuming the isotropic distribution of 
9. A logarithmic divergence near ~ 7r/2 is removed 
by the cutoff \/2d/a. We thus obtain (|1/ cos 6*1)8 = 
(l/7r)ln(2a/d) and 



r = 2 — -^ln(2a/d). 



(C2) 



For evaluation of the winding time Ti we introduce 
a probability q that a particle changes its direction of 
motion along a ring arm after scattering from the ring 
boundary. In the case of diffusion scattering q — 1/2. 
If g < 1/2, the specular reflection prevails. In such a 
situation the time t^^^ = ^^Th plays a role of a mean 
free time for a particle that propagates diffusively along 



a ring arm. The corresponding diffusion coefficient D 
can be evaluated as {IS.x1fj)s/T'fJ^, where (Aa^gy^j)^ = 
{vFTl^^ f{&iTi^e)s is the mean quadratic distance along 
a ring arm that an electron passes during the time t^^^ ■ 
Calculating the average (sin^ 6)s in the same way as 
above we obtain 



D 



Awq 



ln(2a/rf) . 



(C3) 



The distance L passed by a diffusing particle during the 
time r is L = V DT and, for the winding number we 
obtain, accordingly 

{w^), = iV(27ra)2 = DT/{2naf = T/Ti. (C4) 

Finally, we get from this equation and Eq. (C3) 

d 



"'Ms) 



d/ VF^T^{2a/d) 



(C5) 



To find r_L we use definition (61) together with 
Eqs. (52) and (49). A simple algebra gives 



T±=Ti 



47r2 7^(1 + 



(C6) 
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